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Inviscid  modes  of  instability  in  spiral  flow* 

E.T.Y.  Lee**  and  W.  H.  Reid 
Brown  University 
Providence,  Rhode  Island 


Abstract 

The  stability  of  an  inviscid  fluid  between  rotating 
cylinders  with  an  axial  flow  is  considered  under  the  small  gap 
approximation.  When  the  axial  flow  is  small  compared  to  the 
rotational  velocity,  a  new  perturbation  method  differing  in 
several  respects  from  the  conventional  method  has  been  devised. 
It  Is  found  that  the  presence  of  the  axial  flow  has  a  stabiliz¬ 
ing  effect,  and  that  the  correction  to  the  growth  rates  of  the 
pure  rotation  case  is  of  second- order  in  the  axial  to  rotational 
velocity  ratio.  Furthermore,  the  instabilities  are  confined 
to  only  a  finite  range  of  wave  numbers . 
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1 .  Introduction. 

The  effect  of  an  axial  pressure  gradient  on  the 
stability  of  viscous  flow  between  rotating  cylinders  has 
recently  been  studied  both  theoretically  by  Chandrasekhar 
[lr3]  and  Di  Prima  [4]  and  experimentally  by  Donnelly  &  Fultz 
[5]  and  Snyder  [6],  This  problem  is  of  particular  interest 
because  of  the  interaction  it  exhibits  between  the  Tollmien- 
Schlichting  mechanism  of  instability  associated  with  the  axial 
flow  and  the  Taylor-Qortler  mechanism  associated  with  the 
rotational  flow.  As  in  the  usual  stability  theory  for  parallel 
flows  in  which  the  inviscid  form  of  the  Orr-Sommerfeld  equation 
plays  a  central  role,  the  corresponding  inviscid  analysis  for 
the  present  problem  may  be  expected  to  be  of  comparable 
importance . 

In  the  present  discussion  we  will  make  the  usual  small  - 
gap  approximation,  i.e.  we  assume  that  the  spacing  between  the 
cylinders  is  small  compared  to  their  mean  radius.  We  will 
also  restrict  the  discussion  to  velocity  distributions  that 
are  physically  realizable  in  a  viscous  fluid  so  that  the 
angular  and  axial  velocity  distributions  are  given  by 

Q(r)  =  aJi-U-^]  (1) 

and  W(r)  £  6Wmtfl<),  (2) 

where  £  =  (r-R^/d,  ^  and  is  the  mean  velocity 

of  the  axial  flow. 
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For  an  inviscid  fluid,  the  stability  of  a  purely  rota¬ 
tional  flow  is  determined  by  Rayleigh’s  circulation  theorem 
which  states  that  a  necessary  and  sufficient  condition  for 
stability  with  respect  to  rotationally  symmetric  disturbances 
is  that  the  square  of  the  circulation  should  increase  outwards 
from  the  axis  of  rotation.  On  this  criterion,  the  distribu¬ 
tion  of  angular  velocity  given  by  equation  (l)  is  unstable  if 
jt  <  1  and  stable  otherwise .  On  the  other  hand,  the  stability 
of  a  purely  axial  flow  is  governed  by  Rayleigh’s  inflection 
point  theorem  which  states  that  a  necessary  condition  for 
instability  is  that  the  velocity  profile  have  an  inflection 
point  somewhere  in  the  flow  region.  For  the  symmetrical 
velocity  distribution  (2),  this  condition  is  also  sufficient. 

We  would  expect,  therefore,  that  the  stability  or  Instability 
of  a  basic  spiral  flow  would  be  determined  solely  by  its 
rotational  component  as  was  first  suggested  by  Chandrasekhar  [73. 

Since  the  axial  flow  is  Btable  in  the  absence  of  rota¬ 
tion,  we  might  expect  it  to  have  a  stabilizing  effect  on  an 
unstable  rotational  flow  in  the  sense  that  the  growth  rates 
of  the  disturbances,  if  positive,  would  be  decreased  when  an 
axial  flow  is  also  present  and  that  the  instability  would  then 
extend  only  to  a  finite  range  of  wave-numbers .  This  conjecture 
is  borne  out  by  the  present  calculations. 

When  the  axial  flow  is  dominant  compared  to  the  rota¬ 
tional  flow,  however,  the  situation  is  more  complicated  for 
the  governing  equation  then  possesses  a  singularity  in 
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the  region  of  flow.  In  the  total  absence  of  rotation,  the 
roots  of  the  lndlclal  equation  relative  to  this  singular  point 
differ  by  an  Integer  and,  consequently,  one  of  the  solutions 
(commonly  denoted  by  cp2)  has  a  logarithmic  singularity.  When 
rotation  is  present,  however,  the  roots  of  the  lndlclal  equa¬ 
tion  no  longer  differ  by  an  integer  and  the  singular  point  Is 
than  a  regular  one.  Thus,  even  a  small  amount  of  rotation 
must  result  In  a  large  modification  of  the  solution  <p2  In  the 
neighborhood  of  the  critical  layer.  This  aspect  of  the  problem 
will  be  considered  further  in  a  later  paper. 

In  the  existing  non-dimensional  formulations  of  the  full 
viscous  problem  two  parameters  appear:  the  axial  Reynolds 
number 


B  =  V/v 


(3) 


and  the  usual  small-gap  Taylor  number 


4A%  ,4..  wW*d 

- d  £  2(l-ji)(  -  )  jr 


T  =  - 


(4) 


v  -  '  -ll 

In  considering  the  approximate  solution  of  this  problem  by 
asymptotic  methods  it  is  convenient  to  replace  either  R  or  T 
by  a  new  parameter  that  is  independent  of  viscosity.  Hius, 
we  define  the  quantity 

P  s  R//T:  (5) 


When  the  effects  of  rotation  are  dominant,  the  relevant 
parameters  are  clearly  T  and  p ;  conversely,  when  the  axial 
flow  is  dominant  we  would  use  R  and  1/p .  The  related  inviscid 
problems  are  then  obtained  by  formally  allowing  T  or  R  to 
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become  infinite  for  fixed  values  of  p.  Both  of  these  limiting 
procedure  are,  of  course,  singular  since  the  order  of  the 
governing  equations  is  thereby  reduced  from  six  to  two. 

In  all  of  the  existing  work  on  the  viscous  problem  it 
has  been  assumed  that  the  axial  Reynolds  number  1b  small  so 
that  the  effects  of  rotation  are  dominant.  The  present 
inviscid  analysis  will  also  be  confined  to  this  case  where 
p  <<1.  Even  in  this  very  restricted  fora  the  problem  presents 
a  number  of  difficulties  and  it  has  been  necessary  therefore  to 
resort  to  a  perturbation  method. 
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2.  The  governing  equations . 

If  the  steady  basic  flow  given  by  equations  (l)  and  (2) 
is  subjected  to  a  rotationally  symmetric  disturbance  whose  t 
and  z  dependence  is  of  the  form 

exp(pt+ikz),  (6) 

then  the  linearized  equation  for  u,  the  radial  component  of 
the  disturbance  velocity,  can  be  written  In  the  non-dimensional 
form  (cf.  [2]) 

(<j  +iapw)2(D2-a2)u-ia(d  +iapw)(D2w)u  *  -a2«u,  (7) 

where 

D  -  d/d£,  a  -  kd,  and  c  =  p/(-4A£21)^.  (8) 

In  equation  (7),  u  and  w  are  the  non-dimensional  forms  of  the 
angular  and  axial  velocity  distributions : 

«(0  =  l-(l-n)C  and  w  (£)  =  C(l-C).  (9) 

The  parameter  {3,  defined  by  equation  (5),  can  also  be  written 
in  the  form 

p  =  6(wm/a1d)(-4A/D1)“^  (10) 

in  which  there  is  no  longer  any  explicit  reference  to  the 
viscous  problem.  The  boundary  conditions  that  must  be  satisfied 
by  u  are 

u=0  at  £  =0  and  1.  ( 11 ) 

The  purely  rotational  case  (p»0)  has  been  considered  by  Reid  [8] 
who  obtained  an  exact  solution  for  -00  <  <  1  in  terms  of  Airy 

functions . 
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If  we  confine  our  attention  to  the  case  in  which  the 

cylinders  rotate  in  the  same  direction,  i.e.  to  the  range 
0  <  V-  <  then  “(O  lB  of  one  sl8n  611(3 *  to  a  Sood  approxima¬ 
tion,  can  be  replaced  by  its  average  value  um  -  ^(l+y,) .  More 
precisely,  if  we  write 

w(C)  -  wmtl-e(C-|)  ],  where  e  -  2  p-Jj  ,  (12) 

then  it  has  been  found,  in  both  the  viscous  [2]  and  the 
lnviscid  [53  problems  for  the  purely  rotational  case,  that  the 
correction  to  the  eigenvalues  is  only  of  second-order  in  e . 
This  approximation  is,  of  course,  no  longer  valid  when  ji  <  0, 
since  w(C)  will  then  change  sign  in  0  <  C  <  1.  Hie  simplifica¬ 
tion  that  results  from  this  approximation  can  be  readily  seen 
if  we  note  that,  in  a  perturbation  method  based  on  the  small¬ 
ness  of  the  parameter  P ,  the  lowest  approximation  will  be  the 
purely  rotational  case  with  u  a  constant,  the  solutions  of 
which  will  be  simply  sines  and  cosines  rather  than  the  more 
complicated  Airy  functions.  Furthermore,  the  analysis  to  be 
described  can,  in  fact,  be  carried  through  for  the  general 
case  with  but  slight  modifications,  and  it  is  found  that  most 
of  the  important  features  of  the  problem  are  already  exhibited 
in  the  simplified  case. 

Thus,  with  u(£)  replaced  by  its  average  value  u^,  we 
have  the  governing  equation  in  the  final  form 

(d+ia{3w)2(D2-a2)u-ia(c+iapw) (D2v)u  *  -a^)mu.  (13) 

For  p  «  1,  an  approximate  solution  of  this  equation  can  be 
obtained  by  the  perturbation  method  described  in  the  following 
section . 
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3.  The  perturbation  method. 

The  method  to  be  described  here  differs  from  the  usual 
perturbation  methods,  as  described,  for  example,  by  Courant 
and  Hilbert  [9],  In  a  number  of  important  respects  the  most 
important  of  which  is  the  elimination  of  the  crucial  require¬ 
ment  that  the  normal  modes  of  the  unperturbed  problem  be 
complete.  In  the  present  unperturbed  problem,  for  example, 
when  the  cylinders  rotate  in  opposite  directions  (i.e.  jx  <  0) 
the  completeness  question  has  not  been  satisfactorily  answered. 
In  addition  to  the  discrete  spectrum  of  positive  and  negative 
eigenvalues  (with  limit  points  at  +oo)  that  has  been  found,  the 
problem  may  also  possess  a  continuous  spectrum  (though  this 
later  possibility  has  not  yet  been  investigated) . 

To  avoid  this  difficulty,  we  shall  fix  our  attention  on 
a  particular  mode  and  show  that  it  is  possible  to  obtain 
corrections  to  the  eigenvalue  and  eigenfunction  of  that  mode, 
to  arbitary  orders  in  p,  without  any  knowledge  of  the  other 
modes.  Prom  the  point  of  view  of  stability  theory  this  method 
provides  all  of  the  information  that  is  usually  required.  When 
the  cylinders  rotate  in  the  same  direction  and  u  is  replaced 
by  its  average  value  as  in  equation  (13),  the  eigenfunctions 
of  the  unperturbed  problem  are  clearly  complete  and  there 
would  be  no  difficulty  in  applying  one  of  the  usual  perturba¬ 
tion  methods .  Even  in  this  case,  however,  the  present  method 
has  the  important  advantage  that  the  corrections  to  the  eigen¬ 
value  and  eigenfunction  of  a  given  mode  can  be  obtained  in 
terms  of  finite  sums  rather  than  infinite  series . 


562(07)/52 


9 


The  success  of  the  present  method  depends  crucially  on 
being  able  to  solve  explicitly  the  Inhomogeneous  equations 
that  result  from  the  expansion.  This  has,  however,  proved  to 
be  possible  not  only  in  the  present  problem  (for  arbitrary 
rates  of  rotation)  but  also  in  Bisshopp's  work  [10]  on  non- 
rotationally  symmetric  inviscid  modes  in  Couette  flow.  In 
fact,  Bisshopp’s  work  provides  sm  example  in  which  the  "averag¬ 
ing"  approximation  was  not  made  and  the  analysis  was  carried 
through  in  terms  of  Airy  functions . 

Before  we  consider  the  expansion  for  the  solution  u, 
it  is  convenient  to  introduce  the  transformation 

x  »  XjC  where  x1  ■  (1*0 

As  a  result  of  this  transformation  the  value  of  x  correspond¬ 
ing  to  C "1  at  which  one  of  the  boundary  conditions  must  be 
applied  will  depend  on  {3 .  The  function  u  is  now  of  the  form 
u«u(x;a,n,d,p)  with  d=d(a,ji,  p)  but  their  dependence  on  a 
(supposed  fixed)  and  n (which  occurs  only  through  um)  need  not 
be  indicated  explicitly  in  what  follows.  We  now  expand  u  in 
the  form 

u  -  uq(x,cj)  +  Pu1(x,d)  +  p2Ug(x,cj)  +  ...  (15) 

but  do  not  expand  d  until  a  later  stage  in  the  analysis.  On 
substituting  the  expansion  (15)  into  equation  (13),  and  collect¬ 
ing  terms  of  the  same  order  in  p,  we  obtain 
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(D2+1)u0  -  0,  (16) 

(D2+l)un  -  u  2  c.x1,  (17) 

1  0  l-o  1 

2  3 

(D2+l)ug  =  vl^  2  CjX1  +  uQ  Z  djX1,  (18) 

1=0  1*0 

etc.,  where  D  &  d/dx,  c.^  =  ci(a,d,x1),  and  dA  =  d1(a,d,x1). 

To  avoid  going  Into  too  much  detail,  we  will  not  give  the  full 
expressions  for  ci,  d.^  and  the  other  coefficients  that  will 
appear  later;  it  suffices  to  indicate  the  parameters  on  which 
they  depend.  The  boundary  conditions  now  require  that 
Uo+P^i+p  Ug+ . . .  =  0  at  x=0  and  x=x1,  where  x1=x1(d)  and  d,  in 
turn,  depends  on  0.  Thus,  in  contrast  to  the  usual  perturba¬ 
tion  methods,  the  position  of  the  outer  boundary  is  not  fixed 
and  we  cannot  separate  the  boundary  conditions  to  obtain  inde¬ 
pendent  conditions  for  each  u1(x) .  Instead  we  first  construct 
two  linearly  independent  solutions  of  equation  (13)  correct  to 
0(P  ),  and  the  subsequent  expansion  of  d  then  allows  the 
boundary  condition  at  x1  to  be  satisfied  to  0(p2) . 

Two  linearly  independent  solutions  uQi(i=l,2)  of  equa¬ 
tion  (l6)  are  simply  sin  x  and  cos  x.  Due  to  the  linearity  of 
equations  (17)  and  (18),  particular  solutions*  of  the  inhomo¬ 
geneous  equations  can  be  generated  from  the  solutions  of 

(D2+l)fn  =  xnu0  (19) 

■» - 

Hie  solutions  of  the  homogeneous  equation  for  u^  and  u^  need 
not  be  considered. 
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and 

(D2+l)gn  -  x\,  (20) 


which  are,  as  can  be  easily  verified. 


n 


U\  -  §  X1 


n-1, 


u. 


f  1 
2  n-3 


(21) 


and 


g, 


n 


2TrFTJ 


[xn+1u 


_  n(n+l)  n-l  '  +  (n+l)n(n-l)(n-2) 
2n-l  x  uo  2n-I 


(22) 

In  the  above  formulas  n  is  a  non-negative  Integer,  fft  and  gn 
with  negative  subscripts  are  to  be  regarded  as  zero  Identically, 
and  f1  is  to  be  taken  as  xu^/3  instead  of  (xu^/3  -  uQ/6)  since 
(D2+1)u0  =  0. 

Thus,  corresponding  to  uc=uoi,  we  have,  for  a  particular 
solution  of  equation  (17), 


U11  ■  JL  Cnfn*’ 

n=o 

where  fRl  is  the  solution  (21)  of  equation  ( 19 )  when  uQ 
Collecting  terms  we  obtain 


(23) 
is  u0l . 


uli  “  PlWUol  +  P2(x)uol' 


(24) 


in  which 

2  i  3  , 

P2(x)  -  PjjXJ,  P2^x)  "  ^  *  (25) 

and  the  j=Pij(a#d,x1)  are  linear  combinations  of  the  C^'s. 
Equation  (l8)  then  becomes 
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(D2+1)U2  -  u0  +  u; 


4  .5 

I 

J-l  J  v  J« 

where  mj  and  nj  are  combinations  of  c^d^  and  p^. 

particular  solutions  must  therefore  be  of  the  form 

4  5 

Z 

j-1  J  JJ-  3«1 

«j(x)uoi  + 


(26) 


Its 


^1  *  ,?-,mjfjl  +  ^i^Jl 


(27) 


in  which 

6  ,  5  < 

Q^(x)  *  Z  q^jX  »  Qg(x)  **  ^Z^cjgjX  ,  (28) 

and  the  q±j  *  qij(a,c>x1)  are  linear  combinations  of  m^  and  n^. 

We  have  thus  obtained  two  fundamental  solutions  of 
equation  (13) which,  correct  to  0(p2),  can  be  written  in  the 
form 

ui(x)  -  uQl  +  P(Voi+P2uoi)  +  P2(Qluoi+Q2uoi)  +  ‘  ' 

(29) 

in  which  uQl(x)  =  sin  x  and  uq2(x)  ■  cos  x*  Since 

P1(0)  *  Q1(0)  *  0  and  uq1(0)  =  0,  the  characteristic  equation 

is  simply 

ul(xi)  *  0,  (30) 

and  the  required  eigenfunction  is  a  constant  multiple  of 
u^x)  »  sin  x  +  pfP^xJsin  x  +  Pg(x)cos  x} 

+  p  2  {^(xjsin  x  +  Qg(x)cos  x  }  +  ... 


.  (31) 
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We  may  remark  here  that  in  the  general  case  where  w(C) 
is  not  approximated  by  its  average  value,  a  different  trans¬ 
formation  x=x(Cja,^,d )  must  be  used  with  the  result  that  both 
boundary  positions  are  then  functions  of  d(0) .  The  character¬ 
istic  equation  will  also  involve  both  of  the  boundaries,  x 

o 

and  x1  say,  and  the  required  eigenfunction  will  then  be  a 
certain  linear  combination  of  the  corresponding  fundamental 
solutions  u^x)  and  Ug(x),  the  leading  terms  of  which  will  be 
the  Airy  functions  Ai(x)  and  Bi(x)  respectively. 

It  is  at  this  stage  that  we  will  expand  the  eigenvalue 
<J  in  powers  of  p,  and  thereby  find  corrections  to  the  growth 
rates  of  the  unperturbed  problem.  Thus,  we  let 

c(p)  =  do  +pdx  +  p2cj2  +  o(p2),  (32) 

and  we  can  then  obtain  the  corresponding  expansion  for  the 
boundary  position  x^c),  in  the  form 

x1(o)  »  Xx  +  px;[(d0)o;L 

+  P2^xi(^0)o2  +  i  *I<«0>d!}+  (33) 

where  X1  =  x^),  =  (dXj/dd)^,  and  x"(do) 

*  (d2x1/dd2)^_Q.  in  a  similar  manner  we  obtain, 

sin  x^  =  sin  X^  +  px^(dQ)d1  cos  X1 

+  P2{[x^(d0)d2  +  |  x^(d0)<j2]cos  Xj  -  |[x;[(oo)o1]2sin  X1  } 


+  o(p3). 
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cos  x^  “  cos  -  px^(d0)<J^p  sin  X^  +  O(P^), 

dP, (x, )  o 

p.tx,)  -  7,(X.)  +  p [-4--1-]  a,  +  0(f)2), 

p  =0 

and 

Q1(x1)  =  ^L(X1)  +  0(p), 


',(34) 


i 

) 


where  Tr1(X1)  =  [P1(x1)]  and  ^(X)  -  fQ1(x1)]  .  On 

P=*o  A  P=o 

substituting  these  results  into  the  characteristic  equation  (30) 
and  collecting  terms  of  the  same  order  in  p,  we  obtain 

sin  Xx  =  0,  (35) 


and 


(36) 


,  dP0(x, ) 

+  ^l^)*^©)*!  +  J  xi(d0)d2  }•  (37) 


d2  - - r  v 

Xi(cQ) 


Equation  (35)  requires  that  X^  «=  +nrc  (n=l,2,...)  and 
from  equation  (l4)  we  then  obtain 

(38) 

0  1+X 

where 

\  =  a/|X1|  «  a/nn  .  (39) 

Thus,  in  the  absence  of  an  axial  flow,  there  exist  two  modes 
for  each  value  of  n,  one  being  stable  and  the  other  unstable; 
these  modes  correspond  to  the  so-called  "convective  modes"  in 
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the  viscous  treatment  of  the  problem.  By  using  the  expressions 

i .  . 

for  p^,  xi(d0)t  and  <sq,  we  obtain  the  first-order  correction 
to  the  growth  rate  in  the  form 


(40) 


^  (l+x*) 

Since  is  purely  imaginary,  the  correction  to  the  growth 
rate  is  of  second-order  in  p .  Only  one  value  of  d  1  obtains 
for  each  n,  and  it  is  independent  of  the  rotation  rate  p, 
since  um  does  not  appear  in  equation  (40) .  The  required 
expression  for  o  2  can  conveniently  be  written  in  the  form 

*1  d0d2  *  "  IA(1+X2)3  +  B(l+\2)2  +  |]  (4l) 

5X?  4  X?  4 

where  A  =  -|  +  ^-  ^andB»^--^  +  ^5  .  (42) 

Thus,  for  each  value  of  n,  there  are  also  two  values  of  <jg, 
the  signs  of  which  depend  on  those  of  d  ;  d2  also  depends  on 
M-,  but  this  dependence  is  only  through  dQ. 

The  form  given  by  equation  (31)  for  the  eigenfunction 
is  not  e  \tirely  satisfactory  since  the  parameter  p  occurs  not 
only  in  the  coefficients  Pij(a,d,x1),  q1j(a,d,x1),  but  also  in 
the  variable  x=x^  C  •  We  would,  of  course,  prefer  to  have  u^ 
expressed  in  the  form 

U1  "  ul°^  +  ^1^  +  P2ui2^  +  °(P3)*  (43) 

in  which  the  are  independent  of  p.  This  can,  of  course. 
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be  achieved  by  employing  the  same  expansion  procedures  as 
those  described  above  In  obtaining  and  o2 .  The  details  of 
this  calculation  are  somewhat  lengthy  and  in  the  following 
section,  therefore,  we  will  present  only  the  results  for  a 
particular  wave  number,  namely  a«  it . 


562(07)/52 


17 


4.  Results  for  the  first  unstable  mode. 

We  are  mainly  interested  in  the  first  unstable  mode, 
for  which  X1»+x.  Por  this  mode,  the  dependence  of  dQ,  ^  and 
<J2  on  the  wave-number  a  is  shown  in  Pig.  1.  We  note  that  «2 
is  negative  for  all  wave-numbers,  so  that  the  axial  flow  does 
have  a  stabilizing  effect  as  expected.  Over  most  of  the  range 
of  wave-numbers,  the  values  of  d^d.^  and  dg  agree  very  closely 
with  their  asymptotic  behaviors: 


d  =  -I £  a[l  -  +  0(a  j  ],  as  a  0, 


o  x 


(44) 


*  /uZ  [1  -  +  0(a  j  ],  as  a  — >  oo ,  | 

m  2ar  J 

d^  =  [ (l-~ +  ■~|f  +  O(a^)  ],  as  a  — >  0,  1 

\ 

—  -  g—  [(l+^g-)  -  +  0(a  ^)],  as  a  — >  oo  ,1 

xc  a  ' 

x(A_+B.,H)  3A.  +2B,  ,  2  j. 

--  1  1  +  I 


*45) 


as  a  — >  0, 


A..  l,  3Ai  +B-.  ~  2  h 

-  —5^ —  a4[l  +  (— -  |)-»g  +  Ota"1*)], 


HW) 


as  a  — >  cd  , 


where  A1  and  B^  are  the  values  of  A  and  B  when  X-^-x,  namely 


(47) 


Ax  «=  0.03558677  and  Bx  -  0.2193413. 


562(07)/52 


18 


For  this  mode,  we  have  also  obtained  the  corrections 

to  the  unperturbed  eigenfunction,  u^  and  u^,  for  one 

# 

particular  wave-number,  a*rc.  These  results  can  be  written 
In  the  form 

2  3 

-1  /uJJJ  uW  -  /S%  [sin  itC  I  +  cos  z  (^8) 

J“1  J—l 


In  which 


“ll*1  ,  K]_2  *  "1 

^21  *  lt/3  i  Kgg  «  -* 


T 

2*/3  ;  [ 


umup^  *  sin  nC  2^1*2 +  cosuC  , 


in  which 


6)  ,  R12  -  tc2(5^-|),  R^,-  -Tt2^  -4), 


Rl4  *  ~  2),,R15  =  -  ,  Rl6  -  ngp  , 

Rgl  =  -ic(f^  +  2),  Rgg  =  -7t(^2  -6)  ,  Rg3  =  it  -4)  , 


/  (51) 


R  - 

*25  “  ~W  • 


These  results  are  shown  in  Figs.  2  and  3.  We  note  that  u^A' 

(o) 

is  imaginary  while  u£  7  is  real,  and  that  both  are  symmetric 
with  respect  to  £  =  1/2.  Indeed,  this  symmetry  holds  for  any 
mode  X1  »  +nrc  and  any  wave-number  a,  since  both  the  differ¬ 
ential  equation  (13)  and  the  boundary  conditions  are  symmetric. 


This  value  of  a  was  chosen  partly  for  convenience  and  partly 
because  it  closely  approximates  the  critical  value  of  a  in 
the  related  viscous  problem. 
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Such  symmetry  will  of  course  not  obtain  If  m(0  is  not 
approximated  by  Its  average.  [At  this  stage  it  might  seem 
that  the  algebra  Involved  would  have  been  reduced  had  we 
carried  through  the  analysis  on  the  interval  (-£,+£)  rather 
than  (0,1).  This  is  not  so,  however,  for  the  characteristic 
equation  would  not  then  be  as  simple  as  equation  (30)  and  the 
eigenfunction  will  be  a  linear  combination  of  both  u^  and  Ug, 
thus  doubling  the  number  of  polynomials  involved.  Furthermore, 
the  symmetry  condition  provides  a  useful  check  on  the  analysis.] 
From  Figs.  2  and  3  it  may  also  be  noticed  that  the 
functions  iu^1^  and  u|2^,  if  suitably  normalized,  are  very 
nearly  sine  curves.  This  means  that  the  actual  deviation  of 
the  eigenfunction  from  that  of  the  purely  rotational  case  is 
much  slighter  than  might  at  first  have  been  expected. 
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5.  Concluding  Remarks. 

In  the  general  case  In  which  u(C)  is  not  approximated 
by  its  average  value,  two  distinct  types  of  modes  may  be 
present  in  the  purely  rotational  case  [8].  For  0  <  |i  <  1, 
only  the  "convective  modes"  (o2  >  0)  are  present,  whereas  for 
p  <  0,  there  also  exist  "oscillatory  modes"  (d2  <  0) .  Ihe 
same  perturbation  scheme  can  equally  well  be  employed  here, 
though  it  will  be  more  laborious  to  obtain  numerical  results. 
Qualitatively  it  has  been  found  that,  for  the  convective  modes, 
a1  is  imaginary  and  has  two  values,  one  for  each  value  of  dQ. 
For  the  oscillatory  modes,  again  has  two  values,  one  for 
each  but  they  are  real  and  therefore  the  effect  (stabiliz¬ 
ing  or  otherwise)  of  the  axial  flow  on  these  modes  can  be 
determined  without  going  to  second-order  in  P .  The  values 
of  d.^  will,  however,  depend  on  p  in  the  general  case. 

One  of  the  main  qualitative  conclusions  that  can  be 
drawn  from  this  work  is  that,  when  an  axial  flow  is  present, 
there  is  only  a  finite  range  of  wave-numbers  with  positive 

growth  rates.  Indeed,  we  can  easily  show  that  the  equation 
2 

d  +  B  d.  «  0  does  possess  a  positive  root,  a=a^  say,  and 

O  c  O 

hence  for  a  >  aQ  the  flow  is  stable.  From  equations  (38)  and 
(4l),  it  is  more  convenient  to  find  first  the  real  roots  of 

AjO3  +  B^2  -  J  »  0,  (52) 

where  a  ■  1  +  a2A2  and  ^ 


can  then  be  shown 
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that  only  one  root  of  equation  (52)  leads  to  real  positive 
2 

values  of  a!*,  and  from  this  root  we  obtain 
o 


a. 


*3/^  jr  |. 1  -  kJQ  +  p+  °<f >2>  J  •  <53) 


Thus,  the  range  of  wave-numbers  with  positive  growth  rates 
will  be  smaller  the  more  dominant  the  axial  flow  Is  relative 
to  the  rotational  flow. 

The  differential  equation  (13)  possesses  a  singularity 
where  d+iapw  -  0.  Since  a  is  complex,  this  singularity  does 
not  enter  the  flow  region.  Nor  will  the  real  part  of  the 
singularity  lie  in  0  <  C  <  1  if  p  « 1 .  Indeed,  if  we  write 
x  *  C  -  and  denote  the  position  of  the  singularity  by 
x-u^-Xj+ixg  (x1  and  Xg  real),  it  can  be  easily  seen  that,  for 
wave-numbers  a  «  p  _1,  (x*)2  «  (<jQ/2a)p“1  [l  +  0(p)}  ,  while 
for  a  •  0(p_1)  or  larger,  (x*)2  is  at  least  of  the  older  of 
P"2.  Thus  in  any  case  (x^)2 > >  1/4. 

ftie  present  paper  deals  with  only  one  extreme  case  of 
the  problem  of  spiral  flow.  The  other  extreme  is  the  case 
where  p”1  «  1.  As  mentioned  in  the  introduction,  the  presence 
of  even  a  very  slight  amount  of  rotation  will  drastically  alter 
the  nature  of  one  of  the  solutions  (viz.^g)  of  the  case  P_1»0. 
This  problem  may  presumably  also  yield  to  a  similar  perturbation 
method,  but  complexities  of  a  higher  order  must  be  expected. 

Die  general  case  of  arbitrary  p  is  of  course  of  more  Interest 
and  greater  Importance.  A  simple  scheme  such  as  the  one  used 
here  certainly  will  not  apply,  and  a  detailed  study  of  the 
singularities  of  the  differential  equation  must  first  be  made. 
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Table  I 

Uie  values  of  d0,d1#  and  <*2  vs.  a  for  the  first  amplified  mode 

(V*) 


j  a/it 

‘1<’1  i 

/V82 

!  o.i 
i  0.2 

(8:2 
i  0.5 

J 

0.09950  372 
0.19611  61 
0.28734  79 

O.37139  06 

0.44721  36 

-0.03675  954 
-0.07533  730 
-0.11721  79 
-0.16333  94 
-0.21405  29 

-0.00759  0476 
-0.01485  087 
-O.02154  695 
-0.02759  600 
-0.03306  799 

0.6 

;S:S 

0.9 
i  1 

0.51449  58 
0.57346  24 
0.62469  50 
0.66896  47 
0.70710  68 

-O.26922  14 
-O.32838  61 
-O.39092  99 
-O.45620  27 
-0.52359  88 

-0.03814  681 
-0.04307  795 
-0.04812  384 
-0.05353  535 
-0.05953  901 

2 

:  2 
!  6 

0.89442  72 
0.94868  33 
i  0.97014  25 
0.98058  07 

|  0.98639  39 

-1.23818  4 
-1.95276  9 

i  -2.65611  § 

3:8®  i 

;  -0.18907  66 

!  -0.55770  69 

,  -1.38035  0 

I 

coo 

HI 

1  0.99227  79 
0.99503  72 

| 

i  -5.42285  4 
!  -6.79602  4 

1  -16.29982 
-38.81890 
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Table  II 

up-)  and  up)  vs.  C  for  the  first  amplified  mode  (X1  »it) 


!  C 

-vr'tt) 

1 

0 

0 

0 

0.1 

0.442153 

0.721225 

0.2 

0.779177 

1 .26810 

0.3 

0.984531 

1.63120 

0.4 

1.08311 

1.84136 

1  0.5 

1.11072 

1 .91127 

i 


